SCIENTIFIC 

REPORTS 



s 




OPEN 



SUBJECT AREAS: 
FLUID DYNAMICS 
MECHANICAL ENGINEERING 
BIOMECHANICS 

Received 
23 May 2014 



Accepted 
16 July 2014 

Published 
31 July 201 4 



Gray's paradox: A fluid mechanica 
perspective 

Rahul Bale\ Max Hao\ Amneet Pal Singh Bhalla\ Namrata PateP & Neelesh A. Patankar^ 



^ Department of Mechanical Engineering, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA, "^Engineering 
Sciences and Applied Mathematics, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA. 

Nearly eighty years ago, Gray reported that the drag power experienced by a dolphin was larger than the 
estimated muscle power - this is termed as Gray's paradox. We provide a fluid mechanical perspective of this 
paradox. The viewpoint that swimmers necessarily spend muscle energy to overcome drag in the direction of 
swimming needs revision. For example, in undulatory swimming most of the muscle energy is directly 
expended to generate lateral undulations of the body, and the drag power is balanced not by the muscle 
power but by the thrust power. Depending on drag model utilized, the drag power may be greater than 
muscle power without being paradoxical. 
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The motion of a swimming or a flying animal is determined by the interaction between the fluid and the 
animal's body. The efficiency of movement is dependent on the energy expended in this fluid-body inter- 
action. One perspective is that a swimming animal expends muscle power to overcome drag. This perspec- 
tive formed the basis for the analysis on the energetics of dolphin swimming by Gray in 1936, leading to what is 
known as Gray's paradox^ For a long time, Gray's paradox puzzled researchers studying the energetics and 
hydrodynamics of dolphin swimming. Although the paradox is now considered to be resolved^"^, we revisit Gray's 
analysis that led to the paradox. The focus of this work, unlike prior art, is not on the energetics of dolphin 
swimming, rather, it is on questioning the basic premise of Gray's analysis and its generalization to other 
swimming modes. 

Gray's paradox emerged from an attempt by Gray^ to analyze the energetics of two mammalian swimmers. He 
considered a dolphin and a porpoise swimming at 10.01 m/s and 7.6 m/s, respectively. Assuming a turbulent 
boundary layer flow over the surface of the swimmer's body. Gray used an empirical relation for drag that was 
based on a hydrodynamic model for flow over a rigid, flat plate to estimate drag on the swimmer. Muscle power 
data generated during a period of sustained performance by oarsmen was generalized to mammals and used to 
estimate the muscle power generated by the dolphin and the porpoised From this analysis, it was found that the 
power required to overcome the estimated drag was almost seven times greater than the power generated from the 
swimmer's muscles. The imbalance between the power required to overcome drag and the muscle power available 
to swim became the basis of Gray's paradox. 

The approaches to resolve the paradox may be divided into two categories. The first category is based on a 
hypothesis by Gray, himself He hypothesized that dolphins may be employing a mechanism to re-laminarize the 
flow over its body to prevent the flow from becoming turbulent^ Thereafter, much of the work on dolphin 
hydrodynamics and Gray's paradox focused on a mechanism of re-laminarization for drag reduction so that the 
power estimates for swimming were low^'^'^'^. The second category focused on identif)^ing potential errors in the 
data used by Gray. Williams et and Weihs^ argued that measurements of swimming speed, based on on-board 
observations of dolphins, swimming from stern to bow, could be erroneous since the swimming speed of the 
dolphins might have been enhanced by free riding behaviors when swimming close to a ship. In a review on Gray's 
paradox, Fish^ questioned the estimation of power that uses a muscle model based on the sustained performance 
of oarsmen. For steady and sustained work, slow oxidative muscle fibres are used to generate power. Cetaceans, 
which have both fast glycolytic and slow oxidative muscle fibres, may have the ability to use the anaerobic fast 
glycolytic muscle fibres during burst swimming to generate much higher power than that estimated by Gray^ 
There has also been work focusing on the calculation of thrust. Recently, Fish and co-workers have shown that in 
fact, dolphins do produce significant thrust, thus questioning the need to search for drag reducing mechanisms by 
dolphins^. 

In this work we question the basic premise of Gray's paradox on two levels. First, we examine the expectation in 
Gray's analysis that drag power should be less than muscle power. To that end, we present a power balance 
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equation to understand the cascade of power in swimming animals. 
Second, we examine the influence of different drag models on the 
drag power to muscle power ratio. We specifically consider undulat- 
ory swimming as an example problem and inquire whether this ratio 
is greater than or less than unity depending on the drag model used. 
This inquiry goes to the heart of another longstanding question in 
aquatic locomotion - can the force on a swimming body be uniquely 
decomposed into drag and thrust^°? 

Results 

Power budget for swimming. The average power expended by a 
swimming animal over a swimming cycle is 



speed. Here, it is assumed that axial deformation kinematics are not 
significant. Substituting the components of u in Eq. 1 we get^^ 



:|[((7-n)-u]^IA, 



(1) 



where a is the fluid stress acting on the surface of the swimmer's 
body, n is a unit normal on the surface of the swimmer and points 
into the swimmer's body, and the integration is over the surface A of 
the swimmer. The overbar implies an average with respect to time. 
Thus, by definition Pnet is the rate of work done by the swimmer on 
the fluid. 

It is important to note the sign convention (see fig. 1). The velocity 
components are positive when directed in the positive direction of 
the respective axis. The swimming direction is assumed to be in the 
positive x-direction. The lateral oscillation velocity of undulation is 
in the positive and negative y-directions. 

The sign convention for forces needs careful consideration. A 
positive drag force D imposed by the fluid on the swimmer is directed 
in the negative x-direction. Conversely, by action-reaction the swim- 
mer imposes a positive force D on the fluid that is directed in the 
positive x-direction (fig. 1). Similarly, a positive thrust force T 
imposed by the fluid on the swimmer is directed in the positive x- 
direction. Conversely, by action-reaction the swimmer imposes a 
positive force T on the fluid that is directed in the negative x-dir- 
ection (fig. 1). Finally, a positive net axial force Fx = D — T, imposed 
by the fluid on the swimmer, is directed in the negative x-direction, 
and by action- reaction the swimmer imposes a force on the fluid 
which is positive when directed in the positive x-direction (fig. 1). 

In Eq. 1, u is the velocity field accounting for both the velocity of 
body deformations and the swimming velocity. For most undulatory 
(e.g. eels) and oscillatory (e.g. dolphins) swimmers, the velocity of the 
body is predominantly in the swimming direction and the direction 
of lateral body deformations. For the ease of exposition, we present 
equations with only the axial and lateral velocity components. Let u 
and V be the axial and lateral velocity components, respectively. Note 
that the axial component of velocity is the same as the swimming 
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Figure 1 | The sign convention for forces on the fluid by the swimmer, 
and forces on the swimmer by the fluid. Velocity components are positive 
when directed in the positive directions of the respective axis. The velocity 
component u in the positive x-direction is positive, as shown. The lateral 
oscillating velocity of undulation v is in the positive and negative y- 
directions, as shown. 



f = 1 [(^7•n)•wi](iA- 



[{G-n)-v}]dA, 



(2) 



where Piat is the power expended by the swimmer on the fluid in 
deforming the swimmer's body in the lateral direction (defined here 
as the lateral power), and F^ is the axial force from the swimmer on 
the fluid. Note that FxU is the mean power expended due to ante- 
roposterior movement (Pantpos) by the swimmer on the fluid. The 
swimming speed u can be written as the sum of the mean (w) and a 
fluctuation quantity (w'), i.e., u = u + u'. Substitution of this yields 



Pnet = FxU + FxU' + Pia: 



(3) 



Since there is no net external force on a swimming body, 
Fx = F>—f = 0 during steady swimming, where D and f are the 
mean drag and thrust forces, respectively. The drag and thrust on 
a swimmer depends on how it is deflned, and in general, there is no 
unique deflnition^^. This issue will be discussed later. Pnet represents 
the mean power due to the interaction between the fluid and the 
swimmer's body surface. If there are negligible internal losses, all 
the power generated by the muscles would be transmitted to the fluid 
via work done at the surface of the swimmer, i.e., Pmusc^Pnety where 
Pmusc is the mean muscle power. 

Now we recast the power result in terms of the drag power (Du) 
and the thrust power (fu). Since Fx = D—f = 0 and Pmusc^Pmt^ we 
can recast Eq. 3 as 

Pmusc -fu + FxU' + Plat • (4) 

It is clear that Du, — fu, FxU', and Piat are all powers expended by the 
swimmer on the fluid. Since — fu is the power expended by the 
swimmer on the fluid, it follows that fu is the power expended by 
the fluid on the swimmer. Consequently, the following equation is 
obtained 



: + TU^DU + FxU' + Pla, 



(5) 



where the left hand side is the total power gained by the swimmer, 
and the right hand side is the total power expended by the swimmer. 
Eqs. 4 and 5 are two forms of the power budget for swimming. 

Eq. 5 shows that a swimmer gains power not only from the muscles 
(Pmusc) ^ but also gains thrust power from the fluid (fu). This power is 
expended by the swimmer on the fluid by way of drag power (Du), 
anteroposterior fluctuations (FxU'), and lateral deformations (Piat)- 
In fact, since f = D during steady cruising, the thrust power gained 
by the swimmer is expended entirely as drag power, i.e., fu = Du. 
The muscle power is effectively expended in anteroposterior fluctua- 
tions and lateral deformations, i.e., Pmusc = FxU' + Piat- 

Equivalently, the cascade of power may be rationalized as follows. 
The muscle power (Pmusc) is expended directly to deform the body 
(FxU' -\-Piat), which creates thrust power from the fluid on the body 
(fu). The thrust power gained from the fluid is consequently 
expended as drag power from the body to the fluid (Du). 

Application of power budget to swimmers. The power budget will 
be applied to three representative modes of swimming - undulatory 
swimming (e.g., anguilliform swimmer like an eel or lamprey), 
thunniform swimming (e.g., dolphin), and pectoral fin swimming 
(e.g., sunfish). 

Consider an undulatory swimmer (e.g., eel or lamprey) that passes 
a traveling wave along its body or fin in the anteroposterior direction 
to generate movement. In this case, the drag and thrust producing 
regions are intermingled^^'^^. It is well known that animals swimming 
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Figure 2 | Energy budget in undulatory swimming. An undulatory 
swimming animal uses its muscle power (Pmusc) to do work on the fluid 
(Pnet)' Most of this power expended on the fluid is due to lateral 
undulations (Piat)- The power expended on the fluid is primarily 
concentrated in the wake (Pwake)- Finally, all the power expended on the 
fluid is dissipated (viscous) entirely into the fluid during steady swimming. 

at low Reynolds number expend power in producing undulatory 
kinematics Even at finite Reynolds numbers, swimming animals 
expend power to produce lateral swimming kinematics^^'^^'^^. In prior 
work we computed relevant power terms in Eq. 5 for steadily swim- 
ming undulatory animals and found that most of the power is 
expended to generate movement in the lateral direction (Piat), and 
very little or negligible net power is expended to move in the ante- 
roposterior direction, i.e., F^u' was found to be negligible (Eq. 5)^^'^^. 
The power expended to move in the dorsoventral direction was also 
negligible because typically there are no body undulations or move- 
ments in that direction^ Thus, the power balance equation (Eq. 5) 
reduces to Pmusc + Tu^Du-l-Piat. Furthermore, fu = Du and 
Pmusc^Piat during steady undulatory swimming (fig. 2). This was 
found to be the case for low as well as high Reynolds numbers^ \ 

Thunniform swimming (e.g., dolphin) is considered an extreme 
type of undulatory swimming where most of the undulations are by 
the flukes (hydrofoil) at the posterior end of the body. The body 
causes drag whereas the oscillating flukes produce thrust. The drag 
and thrust producing regions are separate with negligible interming- 
ling. The flukes osciflate laterally to create lift forces which are direc- 
ted in the forward direction to obtain thrust^. Prior data from 
carangiform swimmers and analysis based on lift theory implies 
that here too the muscle power is be expended primarily to move the 
flukes laterally (Piat)- The thrust power (fu) is produced at the flukes, 
and the drag power (Du) is expended on the body^. Due to small 
axial velocity fluctuations, F^u' is expected to be negligible com- 
pared to the lateral power expended at the flukes leading to the 
following power balance: Pmusc + Tu^Du + Piat, where Tu = Du and 
Pmusc^Piat during steady swimming. 

In case of swimming caused by pectoral fins (e.g., sunfish), the 
muscle power is expended in moving the pectoral fins primarily in 
the anteroposterior direction^^. Thus, in this case the power balance 
is: Pmusc + Tu'='Du-\-FxU', where fu = Du and Pmusc^Px^^' during 
steady swimming. 

Should drag power be necessarily less than muscle power? Con- 
sider the question whether drag power should be less than muscle 
power, i.e., should Du<Pmusc^ In general, since muscle power is not 
directly expended as drag power, as noted earlier an equivalent 
question to pose is whether thrust power should be less than the 
power to deform the body, i.e., should fu<FxU'-\-Piat^ The latter 
question is preferred because it is the power expended in 
deformations that produces thrust power. The answer to this 
question would depend on how drag and thrust are modeled as is 
explained below. 

Consider scenarios where the drag and thrust regions are phys- 
ically separated on a swimmer's body. For example, the flukes of a 



dolphin are thrust producers whfle the body causes drag. By this 
definition of drag and thrust, the contribution to Du in Eq. 5 is from 
the body region whereas the contribution to fu and FxU' -\-Piat is 
from the flukes^. It can be argued that the thrust power (fu) gained 
back by the caudal fin from the fluid cannot be greater that the rate of 
work done by the fin on the fluid (F^u' + Piat) since there is no power 
source in the fluid. Hence, by this definition of drag and thrust, 

fu <FxU' + Plat, i.e., Du < Pmusc- 

In case of undulatory swimmers where the drag and thrust pro- 
ducing regions are not distinct, the analysis above may not be reas- 
onable. As noted before, there is no clear consensus on how drag and 
thrust should be defined for undulatory swimmers since the drag and 
thrust producing regions are intermingled^°'^^. Hence, there are dif- 
ferent ways of modeling drag. One approach is based on adding 
forward forces from different parts of a swimmer's body and calling 
it thrust. Simflarly backward forces give drag^^'^^. Another approach 
is to consider the total force on the body with respect to time. The 
instances when the net force is forward can be termed thrust and the 
negative portions of the cycle can be termed drag^^'^^. Yet, another 
possibility is to decompose the swimming kinematics into two parts 
and defining one as the drag mode and another as the thrust mode^\ 
Each of these definitions will lead to different models for drag and 
thrust. The utflity of such models is that equating drag and thrust 
models can lead to reduced-order equations to obtain the swimming 
velocity of the animaP\ However, the models could be such that 
Du> Pmusc- An example of such a scenario will discussed in the next 
two subsections. This does not violate any energy principle since the 
corresponding drag and thrust forces are concurrently present on the 
swimmers body^°"^^. Hence, the resultant force is never too large on 
any part of the body. Consequently, the actual power at any point on 
the body never exceeds the muscle power locally. 

The primary issue in generalizing the argument underlying Gray's 
paradox^'^ is now apparent. In general, it is expected that the drag 
power would be lower than the muscle power. Whfle this would be 
true for certain drag and thrust models, like those where drag and 
thrust producing regions are physically separated, this would not be 
true in general for every functional drag and thrust model. If drag 
power is to be considered as the "useful" part of the output power, 
then it should be noted that the total power input to the swimmer is 
Pmusc + f^ (Eq- 5). The implication in Gray's analysis that the muscle 
power is the only source of power gained by a swimmer needs to be 
revised. 

Instead of comparing the drag power to the muscle power 
(Du/Pmusc)^ if one compares drag power to the total power gained 
by a swimmer (Du/ (Pmusc + fu)), then there would be no paradox 
no matter what definition of drag is used. Since drag is always equal 
to thrust {f = D) during steady cruising, Du/ (Pmusc + fu) is always 
less than one. 

Calculation of drag power. In the next subsection we demonstrate 
the effect of different drag models on the drag power to muscle power 
ratio. To do so, one particular case is sufficient. To that end we 
consider an example of undulatory swimming of larval zebrafish. 
This does not imply that our conclusions are specific to only low 
Re swimming like that of larval zebrafish. We choose undulatory 
swimming because there is no unique way of defining drag and 
thrust on an undulatory swimmer^"'^^'^^. 

We consider three different drag models. We use data from 
numerical simulation of live larval zebrafish kinematics and an 
empirical estimate of drag for the three drag models. Numerical 
simulations of steady free-swimming larval zebrafish were carried 
out to compute the hydrodynamic forces on the zebrafish during 
steady swimming. The numerical simulation tool computes the 
swimming speed and forces of fluid-fish interaction using experi- 
mentally extracted zebrafish swimming kinematics as input. The 
three dimensional flow field around the swimming body is fully 
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Figure 3 | The axial force in the swimming direction ( = — on a larval 
zebrafish during steady swimming. The force oscillates about a zero mean. 
Positive values may be interpreted as thrust and negative as drag. Drag and 
thrust according to this definition can be used to compute the drag power 
and the thrust power during steady swimming. 

resolved in these simulations^^. As the experimental kinematics were 
unsteady in nature, kinematics of one cycle of fish's body undulations 
were extracted from the experimental kinematics and repeated over 
time to get steady kinematics. These steady kinematics were used to 
carry out steady swimming simulations. We used experimental data 
for the kinematics of larval zebrafish provided by Melina Hale 
(University of Chicago). The larval zebrafish was 3.9 mm long, 
and the body undulation frequency was 33.3 Hz. The computed 
mean swimming speed was found to be approximately 1 cm/s, which 
was in agreement with experimental data. Details of the numerical 
simulation tool, the zebrafish kinematics, and validation of the zeb- 
rafish simulations can be found in our prior work^^'^°. 

In the first drag model, we estimate the drag using instantaneous 
axial force, computed from the simulation data, on the swimmer. The 
axial force has cyclic variations about zero over a swimming cycle 
(fig. 3). These variations, about zero, can be interpreted as thrust and 
drag (fig. 3). This model results in a drag power of about 0.173 erg/s. 
During steady swimming, we computed that the larval zebrafish 
expends about 1.67 erg/s in the fluid^\ 

For the second drag model, we use an empirical model to estimate 
the drag on a zebrafish similar to the approach used by Gray\ Gray 
assumed that the drag on a swimmer is equal to the towing resistance 
of a rigid body of the same size and shape\ Hence, to estimate the 
drag we replace the zebrafish by a flat plate of same surface area and 
aspect ratio, and move it in stationary water at the wave speed of the 
zebrafish's body undulations (10.3 cm/s). The drag power, for this 
setup, based on an empirical drag model is found to be approximately 
2.84 erg/s, which is an order of magnitude larger than the drag from 
the first model. 

Finally, we calculated drag according a third model proposed by us 
where we have shown that drag and thrust may be separated in a 
meaningful way for undulatory swimmers (by separating the drag 
and thrust producing kinematics) The drag power based on this 
model was found to be about 5.28 erg/s. 

Influence of drag models on the ratio of drag power to muscle 
power. Here we see how different drag models can lead to different 
values of the ratio of drag power to muscle power. Power calculations 
according to three drag models for larval zebrafish were presented 
above. According to the first drag model, we get Du/Pmusc^ 
0.173/1.67 = 0.1. For the second drag model, we obtain 
Dw/Pww5c~2.84/1.67= 1.7. Finally, according to the third model. 



Dw/P^i^5c~5.28/1.67 = 3.16. The results based on second and 
third drag models should be regarded paradoxical according to 
Gray's analysis. However, this is no paradox. The results turn out 
this way simply because of the particular drag model being used. No 
energy conservation principles were violated as discussed in the 
previous subsection. 

Discussion 

A related issue is the efficiency of swimming which is commonly 
quantified by the Froude efficiency^\ It is defined as r]p = 
Du 

^ — = , where P^ake is the mean power in the wake of the swim- 

Tu -\- P-wake 

mer. In undulatory swimming, all the power due to lateral undula- 
tions is dissipated in the fluid, most of which is in the wake of the 
swimmer (fig. 2). Hence, Pwake^Piat^Pmusc- The Froude efficiency is 
meaningful in the context of Eq. 5. On the other hand, in the context 
of the power balance given by Pmusc^Piaty which has no term corres- 
ponding to the so-called "useful" drag power, the Froude efficiency is 
zero^°. Another drawback of using Froude efficiency is the ambiguity 
in the definition of drag and thrust as discussed above^°'^^. For 
example, according to drag models one, two, and three discussed 
earher, the Froude efficiency is f]p = 0.09, 0.63, and 0.76, respectively. 

To summarize, we provide a fundamental fluid mechanical per- 
spective to view Gray's paradox. We note that the muscle energy is 
predominantly dissipated to cause the deformation kinematics of the 
body. The drag power is balanced not by the muscle power but by the 
thrust power. Drag power can be greater than muscle power, depend- 
ing on how it is defined, without being paradoxical. 
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